Astronomy & Astrophysics manuscript no. 4686 


© ESO 2008 


February 5, 2008 





On the distance of GRO J1 655-40* 

C. Foellmii'2, E. Depagne''^ T. H. Dall^ and I. F. Mirabel ^'"^ 

' European Southern Observatory, 3107 Alonso de Cordova, Casilla 19001, Vitacura, Santiago, Chile 
^ Laboratoire d' Astrophysique, Observatoire de Grenoble, BP 53, 38041 Grenoble Cedex 9, France 

^ Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, Campus San Joaquin. Vicuna Mackenna 4860 Casilla 

306, Santiago 22, Chile. 
* On leave from Service d'astrophysique-CEA-Saclay, 91191 Gif-sur-Yvette, France 

Received <date>, Accepted <date> 

ABSTRACT 

Aims. We challenge the accepted distance of 3.2 kpc of GRO J1655-40. We present VLT-UVES spectroscopic observations to estimate the 
absorption toward the source, and determine a maximum distance of GRO J1655-40. 

Methods. We show that the accepted value of 3.2 kpc is taken for granted by many authors. We retrieved in the ESO archive UVES spectra 
taken in April 2004 when GRO J1655-40 was in quiescence to determine the spectral type of the secondary star. For the first time we build a 
flux-calibrated mean (UVES) spectrum of GRO J1655-40 and compare its observed flux to that of five nearby stars of similar spectral types. 
We strengthen our results with the traditional pair method, using published photometric data of the comparison stars. 

Results. We show that the distance of 3.2 kpc is questionable. We determine a spectral type F6IV for the secondary star. We demonstrate in 
details that the distance of GRO J 1655-40 must be smaller than 1.7 kpc. 

Conclusions. The runaway black hole GRO J1655-40 could be associated with the open cluster NGC 6242 which is located at 1.0±0.1 kpc 
from the Sun. Al D < 1.7 kpc the jets are not a superluminal, and GRO J 1655-40 becomes one of the closest known black holes to the Sun. 

Key words, stars: binaries - stars: individual: GRO J1655-40 - distance scale - stars: microquasar - stars: soft X-ray transient - stars: low-mass 
X-ray binary 



1. Introduction 

GRO J1655-40 (a.k.a Nova Sco 94) is a Soft X-ray transient 
(SXT) in our Galaxy (1=345.0°, b=2.2°; R.A.=16''54"'00^ 
Dec.=-39°50"'450- It has been discovered on July 27, 1994 
with BATSE on bo ard the Compton Gamma Ray Observatory 
dZhanget al."l994'). Optical photometry of GRO J1655-40 re- 
vealed a double-wave d modulation with a perio d of 2.6 days 
Jiailvn et alJ Il995bt Ivan der Hooft et all Il997l) . Strong evi- 
dence that the compact object in GRO J1655-40 is a black hole 
was presented by Bailvn et al. ( 1995b) w ho established a spec- 
ti-oscopic period of 2.601 + 0.027 davs. IShahbaz et alJ (Il999l) 
published the masses of the black hole and the secondary star: 
5.5-7.9 Mq and 1.7-3.3 Mq respectively. Given a distance of 
about 3 kpc, the je ts appeared to be supe rluminal: 1.5+0.4 c 
jTingav et al.ll995l) . 1.05 c (.Hiellming & Rupea.l995fl . where 
c is the speed of Ught. Various spectral types of the secondary 
star have been published. We mention i n particular: F3 -F6 IV 
dOrosz & Bailva.1997.) . and F71V-F6III dlsraelian et aljri999.) . 
The distance of GRO J 1655-40 that i s used in the literature is 
3.2+0.2 kpc, although iMkabelet alJ (l2002l) pointed out that a 
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* Based on data obtained at the Very Large Telescope, European 
Southern Observatory, under program ID 073.D-0473(A). 



distance of about 1 kpc cannot be ruled out with the current 
data. The aim of this paper is to show that none of the distance 
estimate methods used for GRO J1655-40 is good enough to 
claim a distance of 3.2 kpc, and that this distance is actually 
based on one very uncertain interpretation. We present unpub- 
lished optical data of GRO J1655-40 to show that the distance 
of the source is certainly smaller than 1.7 kpc. 

2. The problem of the distance 

The distance of GRO J1655-40 usually quoted in the liter- 
ature is the one determined by Hiellming & Ruperj ( 1 19951) 
who present new radio data obtained with the VIA and 
VLBA. The authors build a kinematic model of the radio 
jets of GRO J1655-40 and obtain a distance of 3.2±0.2 kpc. 
The authors do not actually m e asure t he distance d i rectly , 
but refer to iMcKav & KesteveiJ ll 19941) . iTin gav et alJ ll 19951) 
and Har mon et alJ ( 119951) . Their model provides a refine- 
ment of the d istance range obtained b y the cited authors. 
iHar monet alJ lll995l) rely as well on iMcKav & Kestevei] 
( 19941) and lTingav et al.llll995l) for the dista nce. 

The paper by M cKav & Ke steven ( 1994) simply states: "HI 
observations of GRO J1655-40 made with the AT Compact 
Array show solid absorption in the velocity range +10 to -30 
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km s ' , with a further isolated weak feature at -50 km s ' . The 
balance of probabilities is that the distance is around 3.5 kpc, 
unless the -50 km s"' feature is due to an atypical cloud." 
Although providing important information on the absorbing 
material in the direction of GRO J 1655-40, this lAU Circular is 
not a measurement of the distance. 

Independently, Tingav et al. ( 19951 presented new VLBI 
and A rCA data of GRO J 1655-40. Their HI spectrum obtained 
with ArCA (see their Fig. 2) shows a multi-component profile, 
with a weak feature at -50 km s"' too. They intended to in- 
terpret the various components in terms of HI clouds and HII 
regions participating to the mean Galactic rotation, and con- 
cluded that a minimum distance of 3.0 kpc can be inferred. 
They finally give a range of 3 to 5 kpc possible for GRO J1655- 
40. This estimated distance range is relying on the assumption 
that the weak feature observed at -50 km s"' is not moving 
with a peculiar velocity. Although it looks reasonable, this as- 
sumption might simply not be true. These measurements are 
very dependent on the distribution and velocities of various HI 
clouds along the li ne of sight; a d ifficulty that has been claimed 
to be important bv 'Mirabel et alJ (t2002) who note that, in this 
directio n, there are clouds with anomalous velocities up to -50 
km s"' dCrawford et al.ll989l) . i.e. with an amplitude similar to 
that of the weak feature observed in the ATCA spectrum. 

Moreover, Mirabel et al. (2002) have shown that the above 
radio data of iHiellming & Rupenl ( 1 19951) allow to derive with 
certainty only a relativistic upper limit of D < 3.5 kpc, since 
the relativistic time delay of the motion of the ejecta in the sky 
is given by two equations with three unknowns: the angle with 
the line of sight of the jet axis 6, the velocity of the jet /3 - v/c 
and the distance D. Hiellming & Rupen ( 1995) obtained a dis- 
tance of 3.2 kpc by assuming 6 ~ 84°, which actually cor- 
responds only to the upper limit allowed for GRO J1655-40 
by the ab ove equations (see als o Mirabel & Rodriguez 1994). 
However, iMirabel et al ] 1I2OO2I1 claimed that the assumption 
that the jet axis is parallel to the axis of the orbital plane (i.e. to 
the contrary of Hjellming & Rupen) is equally consistent with 
the observations at radio wavelengths. 

A distance of ~3.0 kpc wa s also proposed on the basis 
of optical data bv iBailvn et alJ ( 1995a). The rnethod u sed is 
the "D" method, following Jon ker & NelemansI (l2004 . They 
have first measured the equivalent widths (EW) of Nal-D lines, 
from which they c laim a color exces s E{B - V) of 1.15. They 
used the results of Bar bon et al.l ('1990). who studied the type- 
la supernova 1989B in NGC 3627. These latter authors men- 
tion that the color excess can be properly derived by studying 
the detailed structure of the interstellar Call and Nal-D lines, 
but to do so, a high-resolution and high Signal-to-Noise ra- 
tio spectrum is i ieeded, which is not their case nor the one of 
the spectrum of iBailvn et al. ( 1995a) that has a resolution of 
~ lOA. Consequently, Barbon et al.. Q990 ) determined an em- 
pirical and roughly linear relation between the equivalent width 
of the Nal-D lines and the color excess E(B - V) for six su- 
pernovae. Although not precisely stated, Bailvn et alJ lll995al) 
must have used this empirical relation to determine their value 
of the color excess, and claimed that the results is consistent 
with the EW of other interstellar lines in the optical domain. To 
obtain the distance, they use the standard relations between the 



absorption, the reddening, the colors an d the apparent and ab- 
solute magnitudes (see e.g. lAlleDll973l) . Finally, they conclude 
that the distance of the source is compatible with D ~ 3 kpc, "in 
agreement with the radio observations" of Tingav et al. ( 19^. 

We note that their result is based on the assumption that 
the absorbing material is distributed homogeneously between 
th e source and the obs erver. This has been again challenged 
bv lMirabeletai]ll2002l) who mention in particular that the red- 
dening in this region of the sky occurs in the local arm within 
700 pc from the Sun. 

Similarlv. lBianchini et alJj 19971) measured the color excess 
also by measuring the equivalent widths through gaussian fit- 
ting of the interstellar 5980A Nal-D doublet and the 6613A 
line in their "higher resolution" data. They use the relati onships 
between e quivalent widths and color exce ss given bv lHerbi3 
a975l) and ldella Valle & Duerbeckl (1 19931) to obtain a (mean) 
color excess of 1.13 , in agreement with the value found by 
IBailvn et alJ (l]j)95a). With a resolving power of about 3700, 
not only their spectrum cannot resolve the multi-component 
profile of the NaD lines, but also miss the fact that the lines 
appear saturated (see our Fig.^, and therefore, simply cannot 
be used to measure the extinction. 

On the other hand, Greiner e t alJ ill995t) presented new 
ROSAT X-ray data, from which they infer a distance of "about 
3 kpc". Their method consists of fitting the halo of the ob- 
served radial profile of the source, produced by the scattering 
of the X-rays by the interstellar dust. To fit the radial profile 
of GRO J1655-40 they assume an uniform dust distribution be- 
tween the observer and the source. From this fit, they obtain, a 
value of the effective optical depth at 1 keV of Teff ~ 0.3. 

Fu rthermore, they used the results of IPredehl & SchmittI 
ill 99^) who have studied in details X-ray halos in ROSAT 
sources. From the fractional halo intensity it is possible to de- 
rive the dust column density. The authors show that a good 
correlation exists between the simultaneously measured dust 
and hydrogen column densities, indicating that gas and dust 
must be to a large extent cospatial. Fro m this correlation that 
depends on the optical depth reif(l keVl lGreiner et al] ll 19951) 
obtain for GRO J1655-40 an absorption of Ay = 5.6 mag and a 
hydrogen column density of Nh - 7.0 x 10^' cm"^. However, 
as emphasized by Jonker & Nelemans ( 2004) . their implicit as- 
sumption is that the sight- line for GRO J 1655-40 has the same 
gas-to-dust ratio as the sight-lines for which the relations be- 
tween Ay a nd Nh have been est ablished. To finally compute 
the distan ce. iGreiner et al.l lll99.5t) use the mean extinction law 
given by 'Allen ('1973'): Ay - 1.9 mag kpc"\ As noted by 
Mirabel et al. (2002), GRO J1655-40 is in the Scorpius region 
of the sky which contains rather clumpy optical dark clouds 
in the foreground. Thus the mean extinction law might not be 
accurate enough in this case. 

Similary, lOrosz & Bailvnl (I1997I) could have used the 
"B" m ethod of distance determination Jjonker & NelemansI 
mt they assumed a distance of 3.2 kpc again citing 
Hiellming & Rupen ( 19951), and a color excess of E(B - V) = 
1.3 + 0.1; this latter value being obtained by iHorne et all 



(119961) who used high-quali ty UV spectra obtai ned with HS T. 
However, the cited paper of lHorne etaD(ll996l) is an lAU cir- 
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cular where it is simply stated that "deep 220-nm absorption in 
the HST spectrum suggests E{B - V) = 1.3 mag." 

iHvnes et alJ ( Il998h presented an multi-wavelength dataset, 
obtained with HS T, AAT, RXTE and CGRO. The distance is 
quoted from Hiellmina & Rupen (1995). The authors discuss 
carefully the problem of the extinction. They have measured 
the reddening of the source by using a power-law fit to the UV 
data. As they mention, "no other assumptions about the prop- 
erties of GRO J1655-40 are required, although [they] do need 
to assume an extinction curve." They used the mean Galactic 
extinctiori curve of Seaton ( 1979) and the /?,,-dependent curves 
of jCardelh et al.llll989l) . They claim that both curves give iden- 
tical results, while using the unusual extinction toward cr Sco, 
17° away from GRO J1655-40, gives a much poorer fit, without 
actually mentioning how unusual this extinction is. Their final 
value is E{B - V) - 1.2 + 0.07 mag. However, there is again 
a priori no reason to think that the mean Galactic extinction 
curve is valid in that particular direction. 

The authors also compute d the absolute magnitud e of the 
secondary star They follow lOrosz & Bailvnl ( Il997h for the 
spectral type of the secondary (F5IV). They rescale the mag- 
nitude to the effective radius of GRO J1655-40, and obtain an 
absolute magnitude of My = 0.7 + 0.5, which is very dif- 
ferent fr om the absol ute magnitude of an F5IV star as given 
by, e.g., iGravl (Il992h : 3.2 mag. In particular, it is not clear 
if they assumed that the star was filling its Roche-lobe or 
not. If tr ue, it is rather str a ightforward to compute t he radius 
(see e.g. IPaczvnskilll97ll Llonker & NelemansI Eoo4 . Then, 
iHvn es et al. ' ('1998') corrected the apparent magnitude (17.12) 
of lOrosz & Bailvn (1997) to 17.18 ± 06. A dopting the dis- 
tance of 3.2 kpc bv lHiellming & Rupen '("1995"), and an average 
reddening law (Ry = 3.2), they obtain E(B-V) = 1.25 + 0.17 in 
agreement with previous publications. 

iBeer & Podsiadlowskil (|2002|) reanalyzed th e lightcurve in 
quiescence obtained bv lOrosz & Bailvnl ( Il997l) . In their mod- 
els, the distance is said to be a free parameter, although it is 
not clear if this parameter has been allowed to go as low as 
1.0 kpc. Moreover, in many pla ces they claim that the dis- 
tance of 3.2 kpc of iHielkning & Rupen (.1995) has been used 
"to tighten" their results. They also mention that they obtain 
a reasonable model of the lightcurves with a rather cool disc, 
with a distance "much closer than the distance of 3.2 kpc" (and 
with a higher ;t'^). For that reason (i.e. the distance is too small 
compared to 3.2 kpc), they discarded the model. 

The authors emphasize that their model is very depen- 
dent on the extinction curve. They note that if the extinction 
curve were incorrect, the fitting procedure would compensate 
by choosing diff'erent values for the distance, temperature of 
the secondary star and the color excess E{B - V). Although 
the temperature must be consistent with the spectral type of 
the secondary, the distance and extinction could both be wrong 
and compensate each other The authors claim to obtain values 
for D and E{B - V) in agreement with the (possibly incorrect) 
3.2 kpc distance and (again) a mean Galactic extinction curve. 
However, even reasonable, a such mean extinction curve can- 
not be used to strengthen the confidence of the results. 



3. New observational material 

Our dataset on GRO J1655-40 consists of VLT-UVES spectra 
available in the ESO archive' (prog. ID 073.D-0473(A), RI. 
Rebolo). The spectra were obtained in 2004 with three differ- 
ent central wavelengths: 4025A on the UVES blue side, 5265A 
comprising H/3 and 6300A comprising Ha on the UVES red 
side. The observations were carried out on April 16, 18, 19, 22 
and 25, and also June 18 of 2004. Apart from April 22, two 
spectra of 1440 seconds were obtained for each of the three 
central wavelengths. For every setting, we have retrieved in the 
the archive all possible calibrations: biases, wavelength cali- 
bration frames, flat-fields, order definition frames and format 
check frames. 

All the spectra h ave been reduced u sing the UVES pipeline 
(see e.g. Ballester & Hensbergdl 19951) . The reduction process 
includes bias and inter-order background subtraction (object 
and flat-field), optimal extraction of the object (above sky, re- 
moving sky Unes and rejecting cosmic ray hits), division by a 
flat-field frame extracted with the same weighted profile as the 
object, wavelength calibration, rebinning to a constant wave- 
length step and merging of all overlapping orders. The spectra 
have been shifted to the heliocentric restframe, and smoothed 
with a boxcar of 10. We emphasize that the sky line subtraction 
performed by the pipeline is made through the detection of a 
line that cross the whole slit and is producing a flat contribu- 
tion to the flux across the slit. This important point is discussed 
below. 

For the flux calibration, we have used the method described 
in the UVES quality control pages^. This consists of normalis- 
ing each reduced spectrum by the exposure time, correcting for 
the gain of each CCD and correcting for the atmospheric trans- 
mission^. The final step consists of taking the UVES master 
response curve to provide the F-A flux-calibrated and (atmo- 
spheric) extinction-corrected spectra. We used the latest master 
response curve, obtained in September 28, 2004. It provides 
a flux calibration with a relative error less than 10% (see the 
indicated webpage). 

Since the Signal-to-Noise ratio (S/N) of the blue spectra 
(with a central wavelength of 4025A) was very close to unity, 
they have been discarded from our analysis. The individual red 
"Hy6" and "Ha" spectra have a S/N ratio between 7 and 20 (be- 
fore smoothing). They respectively cover the wavelength range 
from 4785 to 5755A and 5835 to 6805A, both with a disper- 
sion of 0.014740 A/pix. The resolving power of all spectra is 
45 000. This is, in this case, good enough to resolve the Nal-D 
lines, that appear saturated (see Fig.Q). 

According to iBrocksopp et al.1 ( 1200 5V GRO J1655-40 re- 
mained in quiescence since its discovery outburst by BATSE in 
July 1994 until 2005 February 17, when an increase in X-ray 
flux was detected by RXTE. Therefore, we can expect that our 
UVES spectra of GRO J1655-40 were taken when the source 
was in quiescence. We show below that this is indeed the case. 



' http://archive.eso.org' 
^ See http : //www . eso . org/ instruments/uvesTJ^ 



^ See htt p : //www . eso . org/ observin g/df o/qu ality/UVES/qc/response . hi 

for details. 
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Fig. 1. Mean UVES spectrum around the saturated Nal-D lines. 



4. Analysis 

4.1. Spectral separation 

In order to build a clean mean spectrum of the secondary 
star of GRO J1655-40 and determine its spectral type, we fit- 
ted the continuum of individual spectra and attempted a sep- 
aration using the ite rative reconstruction method described in 
iDemers et alJ further enhanced and applied to a Wolf- 

Rayet+O binary by Foellmi et al. (2006). The reader is referred 
to the latter paper for a detailed description of the method. It 
has been in this case applied individually to the serie of spectra 
around H/3 and Ha, and allow us to not only "clean" the spec- 
trum of the source by removing any features unrelated to the 
companion (i.e. a possible contribution of a disk, and/or fea- 
tures that may have been left by the 2D background subtrac- 
tion of the pipeline), but also to build a correct mean spectrum 
by aligning the lines of the moving secondary along the orbital 
motion. 

We started by roughly measuring the positions of the cen- 
ter of the HyS/Ha absorption profile, and transformed them into 
radial-velocities (RVs). We used these RVs to shift the spectra, 
and build a high-S/N absorption-component spectrum. The re- 
sult is a good first-guess spectrum of the stellar object that pro- 
duces these absorption lines (that is, the secondary), whereas all 
unrelated features are smeared out (as much as the RV differ- 
ences are large). It is then shifted back to the original positions 
and subtracted from the original spectra. This provides a first- 
guess of the "second-component" spectrum, clearly revealing 
narrow emission peaks in all individual spectra. 

We used a cross-correlation (CC) technique to measure the 
position of the emission peak that is on top of H/3/Ha, and con- 
tinued the extraction process. We also used CC for measuring 
RVs on the extracted absorption-component spectra during the 
next iterations. We performed in total 3 iterations, since no no- 
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Fig. 2. Results from the spectral separation procedure. On both 
panels, respectively showing the H/3 and the Ha spectra, the 
spectrum of GRO J1655-40 is the lower one, while the ex- 
tracted sky spectrum is shifted up by 0.3 continuum unit. The 
mean S/N ratio in the continuum is about 60 to 70. Important 
lines are labeled. 



ticeable changes were visible on the results of subsequent ones. 
The resulting mean extracted spectra are shown in Fig.|2l 

With the RVs of the absorption lines obtained via the CC 
measurements, we tried to compute an orbital solution with the 
period given by Orosz & Bailvn (1997, P=2. 62157 days). We 
obtained a solution with an orbital K amplitude of 225 + 10 
km s ' in excellent agreement with that of Orosz & Bailyn 
(228.2 + 2.2 km s"'). Our systemic velocity (-88 + 5 km s"') 
is however not in agreement with them (-142.4 +1.6 km s"'). 
Since we have very few points and we did not calibrate in an 
absolute manner our RVs, our value is of course less reliable. 
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The RVs of the narrow emission peaks found in the Ha 
and H/3 extracted spectra show no RV variations at all within 
errors (i.e. < 5 km s"'). Some of these emission lines clearly 
have a multi-component profile. It is very likely that these lines 
are either faint and narrow sky lines missed by the UVES 
pipeline, or formed by some "warm interstellar medium" 
(WIM) (s ee e.g. Revnolds 1985; Domgorgen & Mathis 199^ 
ISembach et alJl200d) . We are confident that they do not be- 
long to the GRO J1655-40 system, because of their constancy 
in RVs, and since their (main) velocity is never consistent with 
the systemic velocity of GRO J1655-40. In particular, we iden- 
tified the lines [O I] /I5577, ^6300, ^6364, [O II] ^6498, ^6560, 
[N II] /16548, /16563, ^6583, /16717 (see Fig.|2}. We also visu- 
ally checked that these lines are actually seen in the 2D images 
of our spectra. Not only do all of them clearly cross the whole 
slit length, but they also appear too faint and noisy to be caught 
by the pipeline. 

We found no evidences in our spectrum of any lines moving 
in anti-phase, nor lines that could not be explained by the sky 
or WIM. In other words, the contribution from a disk around 
the black hole, if any, is not seen in our spectra. 

4.2. The spectral type of the companion star 

In order to study the spectrum of the secondary star in 
GRO J 1655-40, we first performed a spectral synthesis of the 
Ha region. For that we u sed the UVES POP database of spec- 
tra (Ba gnulo et all l2003h ^. taking advantage of the fact that 
this provides spectra taken with the same instrument and cen- 
tral wavelength. We used the slowly rotating star HD 156098 
as a templat e (F6IV, V^5.5 37 mag, Teff=6480 K, logg=3.94, 
Fe/H-0.09: .,Edvardsson et al..l993.). which correspo nds almost 
exactly to the model used bv llsraelian et alJ ( Il999l) . The syn- 
thesis and spectral subtrac tion was done using STARMOD 
llBardetJl985tlMontes et al.t1995. 2000). 

We obtained a very good match to the spectrum of 
GRO J 1655-40, as shown in Fig. |3l A rapid check with other 
F-star templates (see below) confirm that the best match is 
done with HD 15609 8, in excellent agreeni ent with the F6III- 
F7IV spectral types of lShahbaz et alJ ( Il999l) . although we can- 
not exclude an uncertainty of at most one subclass. This un- 
certainty ha^jio_signifl^cant_consequences, as already empha- 
sized by llsraelian et alJ ( 1 19991) in their abundance analysis of 
GRO J1655-40. 

There is no sure indication (within the errors) of emission 
fill-in of Ha (the equivalent width of the "emission" in the 
residual of Ha is 0.1 A). Since we used an averaged spectrum, 
we could expect to see a overall smeared out (over the line) 
contribution from any accretion disk. This is not the case, and 
we can conclude that the disk has a negligible contribution (in 
emission) and that our quiescence spectra are completely dom- 
inated by the secondary. In our best fit, we obtain a rotational 
veloci ty v sin ; = 94 + 8 km s~', in excellent agreement with 
that of llsraeUan eTaP ( 1 19991 v sin / = 93 + 3 km s ' ). 



jhttp : //www ■ sc. eso ■ org/santiago/uvespop7] 
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Fig. 3. UVES spectrum of GRO J1655-40 (dotted line) with 
the best fitting synthesized spectrum using a F6IV stellar tem- 
plate broadened to vsin/ = 94 km s ' (solid line). Below are 
plotted the residuals after subtracting the template from the 
GRO J1655-40 spectrum. The residuals are on the same scale, 
but have been off'set by 0.3 for display purposes. Vertical lines 
denote positions of the Ca lines, and one O fine (shifted down). 

Therefore, we conclude that the secondary star of 
GRO J1655-40 shows a spectral type F6IV, as observed in our 
2004 data. 

4.3. Obtaining the distance from the UVES spectra 

Our approach to determine the distance of GRO J 1655-40 
is based on the comparison of the flux of GRO J1655-40 
in quiescence with that of a nearby single star of the same 
spectral type. This is bas ically the method "A", following 
IJonker & NelemansI (12004'). Since there is an uncertainty on 
this spectral type, we also explore the results with different 
spectral types of the comparison star. We strengthen our results 
by using published photometric data of the comparison stars. 

We have flux-calibrated our UVES spectra of GRO J1655- 
40 following two methods: using the UVES response curve 
and using spectrophotometric standard stars observed during 
the same nights. The comparison between the two methods is 
relevant because the first method is the oretically not produc - 
ing an absolute flux-calibration (see e.g. Bagnulo et al. '20031). 
Since we will compare the flux of GRO J1655-40 with that of 
F stars who have been calibrated with one of the methods only 
(namely the method of the master response curve), we must 
ensure that both methods are equivalent.^ 

The first calibration method uses the master response 
curves, and is described above. The second method was to 

' Technically speaking, it seems not necessary to compare the two 
methods, as long as we use the same master response curve for both 
GRO J 1655-40 and the comparison F stars. However, it gives a much 
stronger confidence to the results if we can show that the flux calibra- 
tion of GRO J1655-40 is actually good. 
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retrieve observations of standard stars with identical setups 
from the archive, taken the same nights of observations of 
GRO J1655-40. These observations were available only for 
the nights of April 16 and 18, 2004, and the spectrophotomet- 
ric standard star used is LTT 3218. To make the comparison 
between the two methods, we used the spectra of these two 
nights only. We might argue that slit losses make the compari- 
son pointless. However, it occurred that the seeing was signif- 
icantly smaller than the slit width (1.6" for the observation of 
GRO J1655-40 and 5.0" for the spectrophotometric standard 
stars) for the two nights of April 16 (seeing ~0.7) and 18 (see- 
ing We are therefore confident that the calibration with 
standard stars is reliable. We finally corrected the spectra for 
the airmass and the (tabulated) atmospheric extinction. Fig. ^ 
shows the flux of the UVES spectrum calibrated with response 
curve, versus the flux of the mean UVES spectrum calibrated 
with standard stars. 

The difference between the two calibrations are roughly 
within one cr for the so-called Ha spectra only, whereas it is 
not the case for the Hfi spectra. When dividing the two cali- 
brated spectra, for each of the H/3 and Ha sets, we obtain the 
respective ratios: 1.10 + 0.01 and 1.00 + 0.02. We note that the 
difference between the two methods observed in the "H/J" re- 
gion is known, as seen in the UVES response curve webpage. 
We thus decided to pursue the analysis with the Ha spectra 
only. 

In order to determine the amount of absorption towards 
GRO J1655-40, we searched for an F6IV star with a known 
distance. For that matter we used again the spectrum of 
HD 156098 that has the crucial advantage of having been ob- 
tained with the same instrument, and central wavelength. The 
spectrum has been flux-calibrated with the UVES master re- 
sponse curve. Moreover, t his sta r has a known HIPPARCOS 
parallax: 19.80±0.72 mas jPerrvm an et al. 19 971) . Since its ap- 
parent magnitude is known (5.537, iTerzan &Bernardll 19811 
IPerrvman et alJl9'97l) . we can compute its absolute magnitude: 
2.0+0.2 mag. 

We rebinned the spectrum of HD 156098 to the ex- 
act dispersion of that of the Ha flux-calibrated spectrum of 
GRO J1655-40. We then divided the two spectra and obtained 
the following mean ratio / = (8.5 + 1.2) x 10""*. We make the 
assumption that the flux from HD 156098 is not affected by 
extinction, which is a reasonable hypothesis given its proxim- 
ity (50±0.2 pc). Therefore, the ratio / of the flux of the F6IV 
star in GRO J1655-40 over that of HD 156098 is diff^erent from 
unity because of only the distance and the extinction toward 
GRO J1655-40. 

We can therefore use the following system of equations: 

mi = -2.5 log(Fi) = Ml -H fl -H 5 log (Di) - 5 (1) 
OT2 = -2.5 log(F2) = M2 + 5 log (D2) - 5 (2) 

where mi and 1112 are the apparent magnitude of 
GRO J1655-40 at quiescence and HD 156098 respectively, q 
is the interstellar extinction towards GRO J1655-40, Mi and 
M2 the respective absolute magnitudes, and Di and D2 the 

^ The values are obtained from the ESO archive DIMM seeing val- 
ues, available at the time of observations. 
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Fig. 4. Comparison between the two calibration methods. The 
flux of the mean UVES spectrum calibrated with standard 
stars is indicated in abscissa, and the flux of the mean UVES 
spectrum calibrated with response curves indicates in ordinate. 
Top panel: spectrum with central wavelength of 5265A {H/3). 
Bottom panel: spectrum with central wavelength of 6300A 
(Ha). The dark line is bissectrix. The gray area shows the 1- 
cr region around the fitted slope through the points. 



respective distances. Fi and F2 indicate the observed flux of 
GRO J1655-40 and HD 156098 respectively, and whose ratio 
is /. The difference of the two equations gives: 

OT2 - mi - 5 log D2 - 5 log Di + M2 - Ml - a (3) 

- Hk) 
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Solving for a we obtain: 



a = 5 log 



D2 1 



^^f) 



+ M2- Ml > 



(5) 



This equation can be used to determine a maximum dis- 
tance of the source. Since we have no precise way to determine 
the absolute magnitude of GRO J1655-40 in quiescence, we 
make the hypothesis that it has the same absolute magnitude 
to that of the comparison star, more or less one magnitude, i.e. 
Ml = M2 + 1 . We note that a does not indicate the standard and 
usual (photometric) absorption in the V band, but rather the ab- 
sorption within the wavelength limits provided by the spectra. 
The above equation is consistent as long as we stay within the 
wavelength limits given by the spectra. These (narrow) limits 
also allow us to assume that the absorption is constant within 
the wavelength range. 

In order to see the dependence of the results with the spec- 
tral type, we made these calculations for five different F stars 
with no known peculiarities (such as a binary or variable sta- 
tus), all retrieved from the UVES POP (and used above for 
spectral synthesis). The five F stars are summarized in Table[2 
where we indicates the star's name, its spectral type, paral- 
lax distance and absolute magnitude. The ratio / between the 
spectrum of GRO J1655-40 and the spectrum of the F star is 
indicated. We are confident that these ratios can be consid- 
ered as constant throughout the wavelength range, since slopes 
of the ratio spectra were extremely small (comprised between 
-1 X 10-^ and2x 10"**). 

The last columns give the ranges of maximum distance of 
GRO J1655-40 computed from the constraint that a > 0, using 
the lower and upper limits of the absolute magnitude of Mi . We 
note that an HIPPARCOS parallax error of 1 mas translates into 
an uncertainty of roughly +3 pc for the distance D2. The source 
of uncertainty on the distance of GRO J1655-40 is therefore 
dominated by that on the flux ratio /. An uncertainty of 20% 
on the flux ratio (taking into account 10% for each of the flux of 
GRO J1655-40 and the comparison star) implies an uncertainty 
of roughly ± 0.2 kpc on the distance Di. 

In Table ^ none of the upper values of the spectroscopic 
ranges of maximum distance exceeds 3.2 kpc. Intere stingly, us- 
ing the absolute magnitude Mi - 0.7 given bv Orosz & Bailvr] 
lll997l) . Eq.|5]impl ies a distance ranging from 0.15 to 1.43 kpc 
for the five F stars. Apart from the range obtained with the star 
HD 16673 (F5V), the overall mean maximum distance, that is 
for obtaining a null absorption (a = 0), is 1.7 kpc (i.e. equiv- 
alent to the distance at null absorption for the star HD 156098 
only, which has the same spectral type as the secondary star in 
GRO J1655-40). We thus consider this value as a strong up- 
per limit to the distance of GRO J1655-40, since the absorp- 
tion is certainly not strictly zero. We emphasize here that the 
only assumption we made is the absence of extinction toward 
the comparison star (but even a small absorption toward the 
sources would not challenge our conclusion), and that the ab- 
solute magnitude of the F6IV star in GRO J1655-40 is that of 
HD156098 within two magnitudes. 

To strengthen these results, we used published photomet- 
ric data of the comparison stars and appUed the pair method. 



More precisely, we can estimate a B magnitude of GRO J1655- 
40 in quiescence of ~ 1 6.65 from the lightcurves published by 
IOros7 & Rail7i]lll997h . Using their V magnitude (17.12), we 
obtain (B - V)cro - 1-53. On the other hand, (B - V) colors 
of the comparison F stars can be retrieved from SIMBAD, and 
assuming that there are unreddened because of their close dis- 
tance (hen ce (B-V)f* - (^-y)F. n), we can use the following 
equations jOlsoJ 1 975t Langl 1 9991 vol II, p33.): 

R^ 3.25 + 0.25 ■ {B - V)f*,o + 0.05 ■ E{B - V) (6) 

Av^R EiB- V) (7) 

where E(B - V) = (B - V)cro - (B - V)f.,o- Using the same 
assumption about the absolute magnitude of the secondary star 
of GRO J1655-40, ranges of maximum distance can be com- 
puted using the standard equation: V - My - Ay - 5 log D-5. 
These ranges are given in the last column of Tabled It can be 
seen that they agree with the ranges given by the spectral flux 
method, although the mean distance (excepting HD 16673) is 
slightly larger: 2.25 kpc. 

5. Discussion 

No astrophysical link has been invoked in the litera- 
ture between GRO J1655-40 and its environment, except 
Mirabel et al. (2002) who mention the possible link of the 
source with the open cluster NGC 6242, since the opposite 
direction of the proper motion vector of GRO J1655-40 that 
they determined with HST is clearly pointing to the cluster. 
This open cluster Ues at a distance of 1 .02 + 0.1 kpc from the 
Sun according to the photometric and radial-velocity measure- 
ments ofXjlushkoya et al. ( 1997). Although the association of 
GRO J 1655-40 with the open cluster cannot be proven with 
our data, the upper limit for the distance of GRO J1655-40 
is fully consistent with the distance of the cluster. Assuming 
th at the link is true, usin g the proper motion values published 
by iMirabel et al .' 1*2005) and the distance between the source 
and the cluster NGC 6242, we can compute the time at which 
GRO J1655-40 would have been ejected from the star clus- 
ter : ~ 6.6 X 10^ years ago. Since the actual age of the clus- 
ter is known: 40.55 x 10^ years (Kharchenko et al. 2005), we 
can deduce that the age of NGC 6242 when the progenitor of 
GRO J1655-40 exploded is roughly 40 millions years. Using 
the canonical isochrones of Pietrinferni et al. (2004) for a solar 
metallicity^, it is possible to obtain the turn-off mass of a clus- 
ter with the given age: between 6.7 and 7.2 Mq (assuming that 
the progenitor of the black hole in GRO J1655-40 has followed 
a single-star evolution). If cor rect, this result favors the lower 
end of the mass range given by Shahbaz et al. '('1999^, and is in- 
compatible with the values obtained by Orosz & Bailvn ( 19^. 
Finally, using the radial veloc ity of the system (-142.4+1.6 
km s"^ lOrosz & BailyDl[l997 t). and assuming the age given 
above, we note that GRO J1655-40 is closer to the Sun than 
NGC 6242 by about 100 pc. 



' jhttp : //www ■ oa-teramo . inaf. it/BASTI/| 
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Table 1. The five F stars used to compute the maximum distance of GRO J1655-40 are summarized. The star's name, its spectral 
type, distance (computed from the HIPPARCOS parallax), absolute magnitude, (B - V) color (obtained from SIMBAD) and the 
ratio of the calibrated spectrum of GRO J1655-40 with that of the star are indicated. The last two columns indicate the ranges of 
maximum distance Di of GRO J1655-40 obtained through the constraint of a > with the two limits of the absolute magnitude 
Ml and for both methods. See text for details. The uncertainty on the distance values is +0.2 kpc. The error on the absolute 
magnitude is obtained from the error on the HIPPARCOS distance. 



Star 


Sp. Type 


Distance 


M 


B-V 




f 


Max. D[ (spec) 


Max. Di (phot) 






(pc) 


(mag) 


(mag) 






(kpc) 


(kpc) 


HD 156098 


F6IV 


50+2 


2.0+0.2 


0.50 


(8.5 + 


1.2) 10-" 


1.08 - 2.71 


1.30-3.27 


HD 16673 


F5V 


21.5+0.5 


4.1+0.1 


0.52 


(9.3 + 


1.3)10-" 


0.46- 1.15 


0.51 - 1.27 


HD 210848 


F7II 


68+4 


1.4+0.3 


0.50 


(11.9 + 


1.7) 10-" 


1.25-3.13 


1.70-4.27 


HD 37495 


F4V 


42+1 


2.2+0.1 


0.50 


(10.2 + 


1.1) 10-" 


0.83 - 2.08 


1.21-3.05 


HD 65925 


F5III 


58+2 


1.4+0.2 


0.40 


(8.4 + 


0.9) 10-" 


1.26-3.16 


1.48 - 3.72 



6. Conclusion 

We have shown that the distance of GRO J1655-40 of 3.2 kpc 
quoted in the literature has been obtained through a refinement 
of a distance range, which was in turn not well established. 
We have determined a spectral type F6IV for the secondary 
star using UVES spectra. By flux-calibrating our spectra of 
GRO J1655-40 during quiescence and comparing them with 
various stars of similar spectral types, we have shown that the 
distance of GRO J1655-40 simply cannot be of 3.2 kpc, and 
is certainly smaller than 1.7 kpc. We consider the possibility 
that the sourc e is associated with the open cluster NGC 6242 
(see Fig. 1 of iMirabel et al]l2002h . which lies at a distance of 
1.0 kpc. With such distance, GRO J1655-40 is not a superlu- 
minal source (the jets speed reaches ~ 0.37c). If the distance is 
confirmed (maybe by the future european satellite Gaia), and 
assuming that the distance of other black-hole binaries are cor- 
rect, GRO J1655-40 would become one of the closest known 
black hole to the Sun. 
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